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In contrast to the inflaton's case, the curvature perturbations due to the curvaton field depend 
strongly on the evolution of the curvaton before its decay. We study in detail the dynamics of 
the curvaton evolution during and after inflation. We consider that the flatness of the curvaton 
potential may be affected by supergravity corrections, which introduce an effective mass proportional 
to the Hubble parameter. We also consider that the curvaton potential may be dominated by a 
quartic or by a non-renormalizable term. We find analytic solutions for the curvaton's evolution 
for all these possibilities. In particular, we show that, in all the above cases, the curvaton's density 
ratio with respect to the background density of the Universe decreases. Therefore, it is necessary 
that the curvaton decays only after its potential becomes dominated by the quadratic term, which 
results in (Hubble damped) sinusoidal oscillations. In the case when a non-renormalizable term 
dominates the potential, we find a possible non-oscillatory attractor solution that threatens to erase 
the curvature perturbation spectrum. Finally, we study the effects of thermal corrections to the 
curvaton's potential and show that, if they ever dominate the effective mass, they lead to premature 
thermalization of the curvaton condensate. To avoid this danger, a stringent bound has to be 
imposed on the coupling of the curvaton to the thermal bath. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Observation of the Cosmic Microwave Background Ra- 
diation (cmbr) anisotropy has now confirmed that struc- 
ture in the Universe originates as density perturbation, 
already present a few Hubble times before cosmological 
scales enter the horizon, with an almost flat spectrum 
[1-4]. The only known explanation for this 'primordial' 
density perturbation is that it originates, during an era 
of almost exponential inflation, from the vacuum fluctua- 
tion of a light scalar field [5] . The primordial density per- 
turbation is predominantly adiabatic in character [2, 14], 
and hence characterized by a single quantity which may 
be taken to be the spatial curvature perturbation £. Ac- 
cording to observation, the spectrum of the curvature 
perturbation (roughly its mean-square) on cosmological 
scales is [1] 

V c = (2 x 1(T 5 ) 2 . (1) 

According to the usual hypothesis, the primordial cur- 
vature perturbation is generated solely by the inflaton 
field, defined as the one whose value determines the end 
of inflation. No other field is involved, and the curva- 
ture perturbation is conserved after inflation ends. As a 
result, the spectrum of the primordial curvature pertur- 
bation in this 'inflaton scenario' is essentially determined 
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by the inflation model alone [15, 16]. The only influence 
of the subsequent cosmology is to determine precisely the 
number N < 60 of e-folds of inflation after our Universe 
leaves the horizon. 

It has been pointed out in Ref. [17] (see also Refs. [18, 
19]) that the inflaton scenario is not the only possible 
way of generating large scale structure from inflation. 
The primordial density perturbation may instead origi- 
nate from the vacuum fluctuation of some 'curvaton' field 
cr, different from the inflaton field. The curvaton scenario 
has received a lot of attention [6, 20-47] because it opens 
up new possibilities both for model-building and for ob- 
servation [48]. 

In the curvaton scenario, the curvaton energy density 
is supposed to be initially negligible so that spacetime 
is practically unperturbed. It becomes significant only 
during some era of radiation domination, when the cur- 
vaton field is undergoing a (Hubble damped) sinusoidal 
oscillation causing its energy density to grow relative to 
the radiation background. (In this paper, we take a to 
be a real scalar field.) While this happens, the curvaton 
energy density p a is proportional to the square of the 
amplitude a of the oscillation, and the fractional pertur- 
bation in both of these quantities is time-independent. 
Finally, the curvaton decays before nucleosynthesis, leav- 
ing behind a conserved curvature perturbation which is 
proportional to the fractional energy density perturba- 
tion that existed before the curvaton decayed, 

1 5p a 2 5a 

C=x r = o r — ■ ( 2 ) 

A p a A a 

In the (quite good [32, 47]) approximation of sudden cur- 
vaton decay, the constant r is the energy density fraction 
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of the curvaton at the decay epoch, 
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(3) 
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The fractional perturbation in a will be related to the 
fractional perturbation of the curvaton field during infla- 
tion by some factor q, 
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where the subscript * denotes the epoch when the cosmo- 

logical scales exit the horizon during inflation. Using the 

i 

well-known expression V% a = H*/2ir, valid for any light 
scalar field, with H the Hubble parameter, one obtains 
the prediction of the curvaton scenario [17, 21] 



i 2 H* 



(5) 



(We are ignoring, in this paper, the possible scale- 
dependence of the curvature perturbation.) 

This prediction depends on the numbers q and r, which 
encode the evolution of the curvaton field between hori- 
zon exit during inflation and the epoch of curvaton decay. 
The evolution is given in terms of the effective potential 
V(a, H) by the well-known equation 



a + 3H a + V = 0. 



and its perturbation 



(5a) + 3H(5a) + V" 5a = , 



(6) 



(7) 



where the prime and the dot denote differentiation with 
respect to a and the cosmic time t respectively. The 
explicit dependence of the potential on H comes from 
the time-dependent fields, which must be present in the 
early Universe in order to generate the energy density. 

The simplest case is the one where the effective cur- 
vaton potential is so flat that the curvaton field actu- 
ally has negligible evolution until the onset of (Hubble 
damped) sinusoidal oscillation. Then 5 = 1, and r may 
be calculated in terms of the curvaton decay rate as de- 
scribed in Ref. [41]. This case may be expected to occur 
[17, 41] if the curvaton is a pseudo-Nambu-Goldstone bo- 
son (pngb), so that it has a potential of the form 



U~ (vm) 2 (l-cos(^)) , 



(8) 



whose flatness is protected by a global symmetry. 

In this paper, we consider the case that the curvaton 
is not a pngb. Instead of Eq. (8), we adopt a form for 
the potential which allows a range of possibilities, 



V(a,H) 



1 



m' 



±cH 2 + g 2 T 2 ) a 



A, 



T n+4 



(n + 4)! M n 



• (9) 



For the curvaton model to work m must be far smaller 
than . The term ±cH 2 represents the effect of those 



fields to which the curvaton has gravitational-strength 
couplings, on the assumption that supergravity is valid. 
It may be positive or negative, and its value will change 
at the end of inflation and again at reheating [52-54]. 
Generically c <~ 1, but during inflation one must have 
c* <C 1 so that the curvaton qualifies as a 'light' field. The 
term g 2 T 2 represents the effective coupling g of the cur- 
vaton to those fields which are in equilibrium at temper- 
ature T after inflation. Finally, we have kept the leading 
term n > in the presumably infinite sum of quartic and 
higher terms. In the case of non-renormalizable terms 
(n > 0), we take M to be the ultra-violet cutoff of the 
effective field theory, say the Planck scale, and assume 
that A„ < 1. 

Using this form for the potential, we study the evo- 
lution of the curvaton field in a very general way. We 
verify that the curvaton contribution to the energy den- 
sity, if initially negligible, remains so until the (Hubble 
damped) sinusoidal oscillation sets in, and we also show 
how to calculate the numbers r and q which enter into 
the prediction for the curvature perturbation. 

The plan of the paper is as follows. In Section II, we lay 
out the stage for our study. In particular, we provide the 
basic equations for the background evolution of the Uni- 
verse and we briefly introduce the curvaton model with 
emphasis on how the curvaton imposes its perturbations 
onto the Universe and on what the requirements for a suc- 
cessful curvaton are. Also, in this section, we discuss and 
justify the choice of the form of the curvaton's effective 
potential. In Section III, we study the evolution of the 
curvaton during inflation and estimate the value of the 
field at the end of inflation as well as the amplitude of its 
perturbations at that time. In Section IV, we investigate 
the evolution of the curvaton field and its density ratio 
with respect to the background density after the end of 
inflation. We find analytic solutions in all cases consid- 
ered, which describe both a(t) and Pa{t). In Section V, 
we employ the results of the previous section to calcu- 
late, in all cases considered, the time when the curvaton 
comes to dominate the Universe or, if it decays earlier 
than this, the corresponding value of the density ratio r 
at the time of decay. We also calculate in detail the value 
of q, relating the curvaton's fractional perturbation at 
decay with its initial fractional perturbation, obtained in 
inflation. In Section VI, we study the effects of tempera- 
ture corrections to the effective potential on the evolution 
of the curvaton and find upper bounds on the curvaton's 
coupling g to the thermal bath, which prevent these ef- 
fects from being destructive. In Section VII, we present 
a concrete model realization of our findings. Assuming 
a realistic set of values for the inflationary parameters, 
we study the model in detail and, by enforcing all the 
constraints and requirements, we identify the parameter 
space for a successful curvaton. Finally, in Section VIII, 
we discuss our results and present our conclusions. 

In this paper, we use units such that c = h = 1 so that 
Newton's gravitational constant is G _1 
rap = 2.4 x 10 18 GeV is the reduced Planck mass. 
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II. THE SET UP 



The curvaton model 



A. The background evolution of the Universe 

The dynamics of the Universe expansion is determined 
by the Friedmann equation, which, for a spatially flat 
Universe, reads 



p = 3{m P H) 2 , 



(10) 



where p is the density of the Universe and H = a/a with 
a being the scale factor of the Universe. 

The dynamics of a real, homogeneous scalar field a 
in the expanding Universe is determined by the Klein- 
Gordon equation of motion given in Eq. (6). It is well- 
known that a homogeneous scalar field can be treated as 
a perfect fluid with density and pressure given, respec- 
tively, by 

Pa = Pkin + V and p a = p kin - V (11) 

with pkin = \o~ 2 being its kinetic energy density. In view 
of this fact, we will model the content of the Universe 
as a collection of perfect fluids with equations of state of 
the form: pi = WiPi, where pi and pi is the density and 
pressure of each fluid respectively and Wi is the corre- 
sponding barotropic parameter. In particular, for matter 
{radiation}, we have w m — {u> 7 = §}• In most of the 
history of the Universe, the total density p is dominated 
by one of the above components. 

The energy conservation (continuity) equation for the 
Universe, reads 



p = -3H{p + p). 



(12) 



The above can be applicable to each of the component 
fluids individually, if only they are independent so that 
their densities pi are conserved, diluted only by the ex- 
pansion of the Universe according to Eq. (12). 

If w is the barotropic parameter of the dominant fluid 
then, for a constant w, Eq. (12) gives [55] 



p oc a 



-3(l+u>) 



(13) 



Using this and Eq. (10), it is easy to find that, for 
w 7^ —1, the Hubble parameter and the density of the 
Universe evolve in time as 



. . 2/ ' 4 / hi \ 

H(t) = —, r and p= — — ■ 

w 3(1 + w) 3(1 + w) 2 V t J 



from which it is evident that a cx t 2 / 3( - 1+w \ 

During inflation the Universe is dominated by the po- 
tential density of another scalar field, the inflaton. Ac- 
cording to this, Eq. (11) suggests that the equation of 
state of the Universe during inflation is w ~ — 1, which, 
in view of Eq. (13), suggests that p ~ const.. Hence (c.f. 
Eq. (10)) the Hubble parameter remains roughly con- 
stant. 



The curvaton is a scalar field a other than the inflaton 
that may be responsible for the curvature perturbation 
in the Universe. One of the merits of considering such 
a field is the liberation of inflation model-building from 
the stringent requirements imposed by the observations 
of the Cosmic Background Explorer (cobe) . Indeed, the 
curvaton model changes the COBE constraint on the en- 
ergy scale of inflation V^/ { 4 into an upper bound [20]. 
Furthermore, in the context of the curvaton model, one 
can achieve a remarkably flat superhorizon spectrum of 
curvature perturbations, which is in agreement with the 
recent data from the Wilkinson Microwave Anisotropy 
Probe (wmap) [1, 2]. Finally, since the curvaton (in con- 
trast to the inflaton) is not related to the Universe dy- 
namics during inflation it can be associated with much 
lower energy scales than the ones corresponding to infla- 
tion and, therefore, may be easily linked with lower en- 
ergy (e.g. TeV) physics. This is why physics beyond the 
standard model provides many candidates for the curva- 
ton, such as sfermions, string axions, or even the radion 
of large extra dimensions. 

According to the curvaton scenario, the curvature per- 
turbations, generated during inflation due to the quan- 
tum fluctuations of the inflaton field, are not the ones 
which cause the observed CMBR anisotropy and seed the 
formation of large scale structure. This is because their 
contribution is rendered negligible compared to the cur- 
vature perturbations introduced by another scalar field, 
the curvaton. In a similar manner as with the inflaton, 
the curvaton field receives an almost scale invariant, su- 
perhorizon spectrum of perturbations during inflation. 
However, at that time, the density of the curvaton field 
is subdominant, and the field lies frozen at some value 
displaced from the minimum of its potential. After the 
end of inflation, the field eventually unfreezes and begins 
oscillating. While doing so, the density fraction of the 
curvaton to the overall energy density increases. As a 
result, before decaying (or being thermalized) the cur- 
vaton may come to dominate (or nearly dominate) the 
Universe, thereby imposing its own curvature perturba- 
tion spectrum. After this, it decays into standard model 
particles, creating the thermal bath of the standard hot 
big bang. 



1. The curvature perturbation 

On a foliage of spacetime corresponding to spatially 
flat hypersurfaces, the curvature perturbation attributed 
to each of the Universe components is given by [21] 



Ci = ~H 



Spt 
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Then the total curvature perturbation £(t), which also 
satisfies Eq. (15), may be calculated as follows. Using 
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that 5p — J^i fyi an d Eq. (12), it is easy to find 

(p+p)c = ^2(Pi+PiKz ■ 



(16) 



Now, since in the curvaton scenario all contributions to 
the curvature perturbation other than the curvaton's are 
negligible, we find that 



Co = C. 



1 



Pa 



1 + W J dec P 



(17) 
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where Co 



2 x 10 5 is the curvature perturbation 



observed by COBE, and C, a is the curvature perturbation 
of the curvaton when the latter decays (or thermalizes). 
From the above, it is evident that (c.f. Eq. (3)) 



Co 



r = — ~ 



P 



The curvaton decays when H ~ T, where 

r = niax{r,j, Tt} 



(18) 



(19) 



with r CT being the decay rate of the curvaton and T T the 
thermalization rate, corresponding to the thermal evap- 
oration of the curvaton condensate. Note that, after the 
curvaton decay, the total curvature perturbation remains 
constant, i.e. Co = Cdcc- 

Suppose that the curvaton, just before decaying (or 
being thermalized) , is oscillating in a potential of the 
form V oc a a , where a is an even, positive number. Then, 
according to Ref. [56], the density of the oscillating field 
scales as 

p a oc aT^fr , (20) 

which means that p = — -^p^pH. Using this, Eq. (15) 
gives 



Co 



a + 2 8a 
6 Y 



(21) 
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where we used that bp a j p a = aSa/a. 

Note also that, using Eq. (12) for the oscillating cur- 
vaton, it is easy to show that 



a + 2 



const. 



(22) 



Thus, in the case of a quadratic potential, the oscillating 
scalar field behaves like pressureless matter, whereas for a 
quartic potential it behaves like radiation. For a (Hubble 
damped) sinusoidal oscillation (a = 2), the above give 
- l(6a/a) dcc (c.f. Eq. (2)). 



2. The curvaton requirements 

Apart from providing the correct curvature perturba- 
tion, a successful curvaton needs to satisfy a number of 
additional requirements. These can be outlined as fol- 
lows: 



• Masslessness: In order for the curvaton field to 
obtain a supcrhorizon spectrum of perturbations 
during inflation, at least on cosmological scales, the 
field needs to be effectively massless, when these 
scales exit the horizon. Thus, at that time, we re- 
quire that 



V" « Hi . 



(23) 



• Misalignment: For the curvature perturbations 
due to the curvaton to be Gaussian, it is necessary 
that, when the cosmological scales exit the horizon 
during inflation, the curvaton field is significantly 
displaced from the minimum of its potential. Thus, 
we require that 



(a — (7 m i n )* > -ff* 



(24) 



• WMAP: According to the recent wmap obser- 
vations, to avoid excessive non-Gaussianity in the 
density perturbation spectrum, the curvaton needs 
to decay (or thermalize) when the density ratio r of 
Eq. (18) is sufficiently large and satisfies the bound 
[3] 



r > 9 x 10" 



(25) 



Obviously, in order for the field not to be thermal- 
ized too soon after the end of inflation, we require 
that its coupling g to the thermal bath, generated 
by the inflaton's decay, is small enough. 

• Nucleosynthesis: Due to the above WMAP bound 
on r, one has to demand that the curvaton should 
have decayed (and not just be thermalized) by the 
time when big bang nucleosynthesis (bbn) takes 
place (at temperature Tbbn ~ 1 MeV), which im- 
poses the constraint 



r„ > h 



BBN 



10" 24 GeV. 



(26) 



C. The effective potential 



The form of the perturbative scalar potential is [57] 



V{a) = ImV + £ 



A,; 



,n+4 



n>0 



(n + 4)! M n ' 



(27) 



where < A„ < 1 and M is some large mass-scale (e.g. 
mp) determining when the non-renormalizable terms be- 
come important. We will consider that the sum is domi- 
nated by the term of the (n + 4)— th order. Wc will also 
take into account temperature corrections due to the pos- 
sible coupling g of the field with an existing thermal bath 
of temperature T. Finally, we will consider corrections 
to the scalar potential due to supergravity effects. These 
corrections introduce an effective mass, whose magnitude 



■5 



is determined by the Hubble parameter H(t) [53]. All in 
all, the form of our scalar potential is 



V(a) = \(m 2 ± cH 2 + q 2 T 2 )o 2 + - K — ^ 
v ; 2 V y ' (n + 4)! M n 



, (28) 



where < c, g < 1 [58] . The masslessness requirement 
demands that, at least during inflation, c, « 1. This may 
be due to some (approximate) symmetry or due to acci- 
dental cancellations of specific Kahler corrections, which 
means that it depends on the particular curvaton model 
[59]. Obviously, the masslessness requirement also de- 
mands m <C if* during inflation. 

Depending on the sign in front of c, the above potential 
may have a global minimum or a local maximum at the 
origin. In the latter case and when y/cH 3> to, gT, the 
potential is characterized by global minima at ±cr min such 
that 



(n + 4)! 



cH 2 M r - 



An 



(cH 2 M n )^ , (29) 



where, in the last equation, we have absorbed A n /(n + 4)! 
into M (for n —* 0, we can set M n — > 4!/A). In this case, 
in order to ensure that the potential remains positive we 
may add a constant Vq > cH 2 a^ in [60]. 

If the local minimum is displaced from the origin then 
one can express the potential in terms of 



o — o — Or, 



(30) 



in which case, for yfcH ^> to, gT, the potential of Eq. (28) 
is recast as 



1 



V(o) = -(n + 2)cH 2 o 2 + J ^ 



x 71+1 

E 

fe=0 



a k o 



n+l-fc 



\ (31) 



where a^ 1 = (n + 1 — k)\ (3 + fc)!. Since the series is fi- 
nite, the above is an exact expression and not an ap- 
proximation. It is evident that all the terms of the 
sum in Eq. (31) are comparable to each-other when 

o — o m in and the sum reduces to u™"^ 4 X)fc=o a k- F° r 
b S> cr m j n , however, the sum is dominated by the term 
of the largest order which is o n+A /(n + 4)!. Obviously, 
when b <C c m i n , only the quasi-quadratic ('quasi' in the 
sense that if = H(t)) term is important. Thus, in all 
cases, the potential can be roughly approximated by 



V(o) ~ cH 2 o 2 + 



T n+4 



M™ ' 



(32) 



where we considered that, when b 3> (T m j n , then b ~ o. 

In the case of +c in Eq. (28), o m [ n — and b = o. 
In the — c case, (T m i n ^ but the complicated structure 
near the local maximum is felt only when the field is 
near b <~ cr m i n • Since we attempt a generic study of the 
behaviour of the field's dynamics, Eq. (32) will suffice 
for y/cH 3> to, gT. It is true, however, that a different 
higher-order term may dominate for different values of 



the field. This is one of the reasons for treating n as 
a free parameter here. In general, we expect that the 
order of the dominant higher-order term increases with 
o. This can be understood as follows. 

Reinstating the A„ coefficients, it is easy to see that 
the ratio of two higher-order terms of different orders is 



Vn K_ fMY 

V n +k A„+fe \ O J 



(33) 



This means that, if all the A„'s are of the same magnitude 
then, for o -c M, the sum of the higher-order terms in 
Eq.(27) will be dominated by the lowest order (quartic 
term), while, for o ^> M, it will be the highest order 
that dominates. Since the sum is infinite, this means that 
V(o 3> M) blows up and the perturbative approximation 
is no longer valid. However, if the A n coefficients are not 
of the same magnitude, then it is possible that the sum 
of the higher-order terms is dominated by a single term, 
higher than the quartic, before blowing up. Increasing o 
may change the order of this term to an even higher one 
for given A„'s as shown by Eq. (33). 



III. DURING INFLATION 

A. Regimes of evolution 

During inflation, our scalar field is expected to roll 
down its potential toward the minimum. This roll may 
be of three distinct types, depending on the initial con- 
ditions of the field or, more precisely, on where exactly 
the field lies originally in its potential. Below, we briefly 
discuss these regimes assuming that, for all practical pur- 
poses, the Hubble parameter during inflation is almost 
constant. We also assume, for simplicity, that the mini- 
mum of the potential lies at the origin (if not one should 
simply set o — > b) and consider the potential of Eq. (32). 
Without loss of generality we take o > 0. 

• The fast-roll regime: The field is in this regime 
when it lies at a region of its potential such that 
V" » H 2 . Due to the masslessness requirement 
during inflation, we expect c, <1. Thus, the field 
may be in the fast-roll regime only if it lies into 
the higher-order part of its potential. Then it is 
easy to see that the border of the fast-roll regime 
corresponds to V"(o IT ) ~ if* , where 



+2 ^ Ct ™+v; 



(34) 



o b ~ (H 2 M n )- 

with (T^ in = (Tmin (if*)- 

Once in the fast-roll regime (i.e. for o > o IT ), the 
field rapidly rolls until it reaches o IT . After this, the 
motion of the field is no longer undcrdamped but 
feels instead the drag of the Universe expansion. 
Any kinetic energy assumed during the fast-roll will 
be depicted very fast after crossing 0{ T as can be 
easily seen from Eq. (6) [61]. 
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• The slow-roll regime: This regime corresponds 
to V" <C H 2 . The field's motion is overdamped by 
the Universe expansion and Eq. (6) can be approx- 
imated as 



3H& 



-V'(o). 



(35) 



which is the familiar form of the inflaton's equation 
of motion in slow-roll inflation. 

• The quantum regime: Classically, the field will 
continue to slow-roll until it reaches the minimum 
of its potential. However, an effectively massless 
scalar field during inflation undergoes particle pro- 
duction by generating superhorizon perturbations 
due to its quantum fluctuations. In that respect, 
the quantum fluctuations may affect the classical 
motion of the field if the generation of the pertur- 
bations becomes comparable to the slow-roll motion 
of the field. Let us elaborate a bit on this issue. 

In a region of fixed size if" 1 (horizon), the effect 
of the quantum fluctuations can be represented [62] 
as a random walk with step given by the Hawking 
temperature So = Th = H*/2ir per Hubble time. 
The slow-roll motion is comparable to this "quan- 
tum kick" when 5o <~ o/H*, or equivalently when 
o — oq, where a is given by Eq. (35) and oq is 
implicitly defined by the condition 



(36) 



(Here a corresponds to region of fixed size H~ x ). 
One can understand the above also as follows. The 
energy density corresponding to the quantum fluc- 
tuations, which exit the horizon, is p qm ~ Tfj <~ 
[63] . The quantum regime begins when this energy 
density is comparable to the kinetic density of the 
slow- roll pkin = \o~ 2 . Comparing the two, one finds 
the same condition, Eq. (36), for the border of the 
quantum regime oq [64] . 

After the onset of the quantum regime, the coherent 
motion of cr* ceases. Instead, the field configura- 
tion spreads toward the origin by random walk. As 
a result, the mean value at horizon exit of the cur- 
vaton field in our Universe (denoted by o„) may lie 
anywhere within the range |<r*| < oq, which means 
that, typically, its value would be ~ oq. The 
same is true if the field finds itself already into the 
quantum regime at the onset of inflation. 

Thus, from the above, we can sketch the typical evo- 
lution of the field during inflation. Since the curvaton is 
expected to have negligible energy density when the cos- 
mological scales exit the horizon during inflation, if we 
demand that initially p a ~ Vj n f , this means that the field 
must originally lie in the steep (and curved) part of its 
potential (where it cannot act as a curvaton) and, there- 
fore, it will engage into fast-roll. Very soon, however, 
fast-roll sends the field below <7f r and the curvaton enters 



its slow-roll regime. Whether it remains there until the 
end of inflation or not depends on the duration of the in- 
flationary period. If inflation lasts long enough, the field 
will slow-roll down to oq, where it will be "stabilized" 
by the action of its quantum fluctuations. The quantum 
drift may further reduce <r* somewhat but, since in the 
quantum regime the field is oblivious of the potential V, 
there is no real motivation to drag it towards the origin. 
Still, er„ <C oq can be selected anthropically, i.e. by as- 
suming that we happen to live at a special place in the 
Universe. 



B. Initial conditions at the end of inflation 

The above enable us to speculate on the initial condi- 
tions of the curvaton after the end of inflation. Indeed, 
given sufficiently long inflation, we expect cr on d ~ oq and 
CTond ~ because the coherent motion of the field ceases 
once it enters into the quantum regime. However, if in- 
flation ends before the field reaches the quantum regime, 
we have cr cn d > oq , where the subscript 'end' denotes the 
end of inflation. Thus, we see that 



(37) 



Using Eq. (36) and also that V ~ oV" for a pertur- 
bative potential, it is easy to see that 



H 2 s 

r^jj ) H * > H * . 



(38) 



Since, in all cases, o > er on( j > oq, we see that the mis- 
alignment requirement for the curvaton is naturally guar- 
anteed by masslessness. 

Let us estimate oq using Eq. (32). Whether oq lies in 
the quasi-quadratic part of the potential or not depends 
on how large is the value of . Indeed it is easy to see 
that 

/ #*/c* c < c* < I 

° Q ~ 1 H*/c ~ (H*M n )^ c, < c ' ' ' ; 



where 



c Q = (H*/M)^. 



(40) 



From the above, since Co, <C I during inflation, it is ev- 
ident that the misalignment requirement is, quite gener- 
ically, satisfied. Using Eqs. (34), (39) and (40), it is easy 
to show that, when c* < Co, 



— ~ < I 



(41) 



and, therefore, fast-roll always ends before reaching the 
quantum regime. 

We should point out here that c* during inflation can- 
not be arbitrarily small. Indeed, due to the coupling 
g between the curvaton and other fields, one expects a 
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contribution of the form g 2 Tfj to the effective mass dur- 
ing inflation, which introduces the bound c* > g 2 . Thus, 
during inflation, 



g 2 < c* < f . 



(42) 



Note, also, that, if c* < (m/ii*) 2 , the soft mass m dom- 
inates the quasi-quadratic term of the potential. 

In the — c case, the above are valid with the substi- 
tution a — > a. If c* < c then ctq ~ <5q 3> cr min . How- 
ever, if Co < c* <C 1 then erg <C <7 m in, which means that 
<tq <~ cr m i n in this case. 



C. The amplitude of perturbations 



is the decay rate of the inflaton field and the subscript 
'reh' corresponds to reheating. The temperature at this 
moment is the so-called reheat temperature T ro h. After 
reheating, the Universe enters a radiation dominated pe- 
riod. 

In the following, we will use Eq. (6) to follow the evo- 
lution of the curvaton field after the end of inflation. At 
the first stage, we will assume that ^fcH ^> m, gT so that 
the potential may be approximated by Eq. (32). Because, 
after the end of inflation, there are no perturbations gen- 
erated due to quantum fluctuations to inhibit the field's 
motion, the curvaton evolves entirely classically, follow- 
ing Eq. (6). Its evolution depends on which of the two 
terms in Eq. (32) is the dominant. Below we study both 
cases individually. 



We can now estimate the amplitude of the curva- 
ton's almost scale invariant supcrhorizon spectrum of 
perturbations. The evolution of 8a after horizon cross- 
ing is determined by Eq. (7), which, for V" <C if* and 
< (6a) < if 2 , gives 



A. The quasi-quadratic case 

Suppose that it is the quasi-quadratic term that dom- 
inates the potential in Eq. (32). Then we can write 



5(j ~ Sao exp 



I / V" 



3 Vii 2 



if* At 



ii*e 



-fjAN 



(43) 



where Sao = H*/2ir at horizon exit, fj = V"/3H 2 and 
AN = if* At is the elapsing e- foldings of inflation. For 
the cosmological scales, AN < 60 by the end of inflation. 
Thus, for V" /Hi < 0.05, we see that da w 5a ~ if*, i.e. 
the perturbation remains frozen. Thus, in view also of 
Eq. (39), for the amplitude of the perturbation spectrum 
we find 



Sa 
a 



if* 



if* 



< — ~ max{c*, c } . (44) 

cnd CTcnd a Q 



V{a)^-cH 2 (t)a 2 , 



(46) 



where, for simplicity, we have dropped the hat on 
a. Inserting the above into Eq. (6) and substituting 
dt = (1 + w)tdr, we obtain 



d 2 a , da 4c 

^ + (1 - W) d^ + J^ - 



(47) 



The above admits oscillating solutions if c is larger that 



c x , where 



V^= -(l-w). 



(48) 



Here we should mention the fact that, while fj 7^ and 
if* 7^ during inflation leave the amplitude of the per- 
turbation spectrum largely unaffected, they do affect the 
slope of the spectrum, causing departure from scale in- 
variance. Indeed, in Ref. [21] it is shown that the spectral 
index of the curvature perturbation spectrum is given by 



n s = \ + 2{f)-e), 



(45) 



where e = -(if/if 2 )* < 1. 



IV. AFTER THE END OF INFLATION 

After the end of inflation, the inflaton begins oscillat- 
ing around its vacuum expectation value. These coherent 
oscillations correspond to massive particles (inflatons), 
whose energy density dominates the Universe. Thus, af- 
ter inflation ends the Universe enters a matter dominated 
period. The inflatons eventually decay creating a thermal 
bath of temperature T. The Universe becomes radiation 
dominated only after the decay of the inflaton field (re- 
heating), which takes place at if rc h ~ ^ini, where Ti n f 



1. Case c > c x 

Assuming negligible initial kinetic density and after 
some algebra, it can be shown that, in this case, the 
solution for a{t) is 



cos 



VcZiEin (t_ 
3(1 + w) U \t 



3(1 -w) . 
—^-=== sin 

4Vc - c x 



2 f^lnfl 
3(1 + to) V*o 



+ 



(49) 



It is evident from the above that, when c — > c x , the os- 
cillation frequency reduces to zero. 

Since V cx H 2 a 2 cx (a/t) 2 and p kin cx a 2 cx (a/t) 2 and 
in view of Eq. (14), we find that 



Pa cx a -K 3 +-). 



Comparing the above with Eq. (13), we obtain 



^cxa-i^, 



(50) 



(51) 
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FIG. 1: Illustration of the scaling solution in the quasi- 
quadratic case, where m e ft ~ \/cH(t). For ti < ts < i3, the 
curvaton field a rolls down a potential which changes in time 
so that the field's value gradually reduces towards er m in. 



which means that the density ratio of the curvaton is 
decreasing and cannot lead to curvaton domination. 

Note here that, although during inflation masslessness 
requires c* <C 1, this is not necessarily so after inflation, 
when c may assume a different value. In fact, unless 
the (approximate) symmetry that enables masslessness 
persists after the end of inflation, we would expect c ~ 1 
and, therefore, c > c x is quite plausible. 



2. Case c < c x 

In this case, Eq. (47) does not admit oscillating so- 
lutions. Instead, for c < c x and negligible initial kinetic 
density, one finds the scaling solution (see Fig. 1) 



a(t) = cr 




(52) 



The second term in the square brackets soon be- 
comes negligible. Again V cx pkin oc (a/t) 2 , which means 
that [65] 



-|(3+u))+2 v 'cx-c 



(53) 



From the above, it is evident that, when c — ► c x the 
scaling of p reduces to the one given by Eq. (50). Using 
Eq. (13), we now find for the curvaton's density ratio 



The above again shows that the curvaton's density ratio 
is decreasing and cannot lead to curvaton domination. 
To make this clearer, consider that, for c — > 0, Eq. (48) 
suggests that p a /p —> const., which means that, for all 
c > 0, the exponent in Eq. (54) is negative. 

The solution for c = c x , with negligible initial kinetic 
density, is easily found to be 



a{t) = ctq 



1 



1 



hi 



Ignoring the logarithmic time dependence, the above re- 
sults in the same scaling for p a as in Eq. (50). 

From Eqs. (49) and (52), we see that 6a/a = const, 
whether the field is oscillating or not. In that sense was 
Eq. (52) called a 'scaling' solution. Thus, in view also of 
Eq. (21), in the quasi-quadratic case we find 



6a 
a 



6a 
a 



(56) 



B. The higher— order case 

This case corresponds to a ~ a 3> a mnl , when the 
scalar potential in Eq. (32) can be approximated as 



V{a) 



,n+4 



(57) 



In contrast to the previous case, where the effective 
mass was determined by the Hubble parameter, in this 
case the roll of the field does not necessarily commence 
immediately after the end of inflation. Indeed, since at 
least during the latest stage of inflation (when the cos- 
mological scales exited the horizon and afterwards) the 
masslessness requirement demands V" <C H, just after 
the end of inflation, the field originally lies in the poten- 
tial at a point where its motion is largely overdamped 
[66]. Consequently, the field remains frozen at the value 
"end- One can easily show this as follows. Considering a 
negligible initial kinetic density, the motion of the field 
is determined by Eq. (35), which gives 



n+2 



1 



1 



9(1 +w) \n + 3 



1 - 



H 2 



YL 

H 2 



(58) 



Thus, as long as Vq H 2 , we see that a w ao, i.e. the 
field remains frozen. However, while H(t) decreases there 
will be a moment when V"(cr cn d) ~ H after which the 
field unfreezes and begins to roll again down its potential. 
Note that, as expected, a similar result applies to the 
perturbation of the field 6a. For example, solving Eq. (7) 
for the matter era following the end of inflation, one finds 



5o = 5ao (£) cos 



3 Ha 
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which, in the regime i 
5a * 5* [1 + 6*0. 



becomes 

Thus, for the amplitude of 
the perturbation spectrum we expect 



6a 
a 



6a 
a 



(60) 



end 



where the subscript 'osc' denotes the onset of the a os- 
cillations. 

As demonstrated in Ref. [56], the energy density of 
an oscillating scalar field in an expanding Universe with 
potential V cx a n+4 scales as 



Pa oc 



(61) 



However, when the scalar field unfreezes, it does not al- 
ways engage into oscillations. Indeed, it is easy to show 
that Eq. (6) with a scalar potential of the form of Eq. (57) 
has the following exact solution: 



Csolu(0 : 



2M" 



(l-tu)(n + 4) -4 
(n + 4)(n + 2) 2 (1 + w) 



(62) 



Of course, in order for the system to follow this solu- 
tion, it has to have the correct initial conditions. How- 
ever, in Ref. [67], a stability analysis has shown that 
the above solution is an attractor to the system for any 
n > n c , where 



1 + 3w 
_ 



w 



(63) 



Attractor solutions are lethal for the curvaton scenario. 
Indeed, one of the principle characteristics of attractor 
solutions is that they are insensitive to initial conditions. 
In other words, if the system is to follow an attractor 
solution, all the memory of the initial conditions will be 
eventually lost. This is a disaster for the curvaton be- 
cause it means that the superhorizon spectrum of per- 
turbations will be erased. To understand this better, 
consider the above solution cr solu (t) . If this is an attrac- 
tor solution then, after the field unfreezes, it will tend to 
approach this solution, i.e. a(t) — ► a so \ u (t). However, if 
we consider a perturbed value of the field a = a + Sa then 
again a(t) — ► cr so i u (£), which means that 6a — > 0. Thus, 
the perturbation spectrum (6a/a) en d, present at the end 
of inflation, will rapidly diminish so that (<5cr/cr)dcc — * 
and, therefore, from Eq. (21), we see that ( a — > 0. 

Consequently, in order to avoid the above catastrophe, 
we have to enforce the attractor constraint : 



n < n c , 



(64) 



i.e. we forbid the curvaton field to lie at the end of infla- 
tion deep enough in the higher-order part of its potential 
so that the attractor constraint is violated. 

The attractor constraint is more stringent in the 
matter dominated era (before reheating) because then 
n c = 2. Thus, if the curvaton unfreezes before reheating 



is completed then it should not lie into the higher-order 
part of V beyond quartic terms. Similarly, for the radia- 
tion era, n c — 6 so the constraint is relaxed up to terms 
of order 8, if the field is to unfreeze after reheating. 

Of course, one can imagine that, even if the initial value 
after the end of inflation does lie in the forbidden regime, 
after unfreezing there is a chance that, while approaching 
the attractor, the field manages to roll down to parts of 
its potential where the attractor constraint is not violated 
before its perturbation spectrum becomes entirely erased. 
Although this is, in principle, marginally possible it does 
not result in anything more than effectively making the 
borders of the forbidden zone somewhat "fuzzy" . Thus, 
in the following, we will disregard this possibility and 
consider the attractor constraint as strict. 

The attractor constraint has a number of consequences 
on the evolution of the field. For example, from Eqs. (57) 
and (61), one finds for the amplitude a of the oscillations 



a oc a 



6 

n + 6 



n + 6 V n c 



+6 
-2 



a±6\ 



(65) 



where we used also Eqs. (29) and (63). The above shows 
that the attractor constraint in Eq. (64) guarantees that 
the ratio a/cr m in is increasing and, therefore, the oscil- 
lation amplitude will never drop as low as cr m i n . This 
means that once the higher-order term is dominant it 
will continue to be so without ever the quasi-quadratic 
term in Eq. (32) taking over. One can also show this as 
follows. The effective mass for the potential of Eq. (57) is 
TOg ff = V" <~ a n+2 /M n . Then, using the above, one can 
show that 



TOeff OC a 



TO off OC 



HimVi^e), (66) 



from which it is evident that, due to the attractor con- 
straint, ml s decreases slower than cH 2 and, therefore, 
the quasi-quadratic term in Eq. (32) will never take over 
once the field finds itself in the higher-order part of its 
potential at the end of inflation. 

From Eqs. (13) and (61), we find that the density ratio 
of the oscillating curvaton scales as 



^oca- 6 (^)+3( 1 +-). 
P 



(67) 



It is easy to see that the exponent in the above is non- 
positive for n > — 2( 1 1 Z 3 Jf ); which is valid both in the 
matter and the radiation eras for all n > 0. Thus we 
conclude that the oscillating curvaton in the higher-order 
case is unable to dominate the Universe, as in the quasi- 
quadratic case. In fact, we have found that the scaling of 
the curvaton's density is such that, for < n < n c , it is 
decreasing slower than the quasi-quadratic case but not 
as slow as the background density of the Universe. 
It is easy to show that the solution in Eq. (62) scales 

as 



pAa solu )oca-^(m). 



(68) 



Comparing this with Eq. (61), one can see that the scal- 
ing coincides in the border case, when n — n c . As shown 
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in Ref. [67], in this border case, even though the energy 
density of the scalar field does scale according to the at- 
tractor solution, the field engages into oscillations. How- 
ever, as we show in Sec. VD, its perturbation spectrum 
is not preserved and so n — n c cannot be allowed. 

For the oscillating curvaton, we have 5a j a = const., 
since both a and a + 5a undergo the same scaling given 
by Eq. (65). This fact, in view also of Eqs. (21) and (60), 
suggests that 



5a 
a 



(69) 



end 



which is the same result as in the quasi-quadratic case 
(c.f. Eq. (56)). We elaborate more on this point in 
Sec. VD. 

Therefore, in both the quasi-quadratic and in the 
higher-order cases, we have shown that there is no sig- 
nificant damping of the curvaton's curvature perturba- 
tion after the end of inflation [68] . Consequently, using 
Eq. (44) , we find the range for the amplitude of the cur- 
vature perturbations of a: 



Co < Ca < max{c»,c } . 



(70) 



We have also shown that, in both the quasi-quadratic 
and in the higher-order cases, the curvaton is unable to 
dominate or even to increase its density fraction with 
respect to the background density of the Universe. This 
means that for the curvaton scenario to work the soft 
mass m should take over before the curvaton decays (or 
thcrmalizes). In Sec. V, wc investigate this further. 



C. The — c case revisited 

Before concluding this section some reference to the 
— c case in Eq. (28) is called upon. Due to Eq. (65), a 
negative sign in front of c has no effect on the dynamics of 
the higher-order case, because the quasi-quadratic term 
remains always subdominant. In contrast, it may have 
an important effect in the quasi-quadratic case. 

The minus in front of c may exist in inflation and turn 
into a plus after inflation ends. This has no other effect 
than setting the value of cr en d to 



&end = fend + C n 



end 



(71) 



where cr™^ = cr m i n (i? en d). The above is so because, in 
the quasi-quadratic case, <7 cn d < cr m i n . Another possibil- 
ity is that the minus existed there both before and after 
the end of inflation, or that it appeared after the end of 
inflation. Again Eq. (71) applies [69]. 

However, now that the quasi-quadratic term, after the 
end of inflation, has a negative sign, the evolution of the 
field may be affected. Let us study this issue. 

In this case, Eq. (46) becomes 



V(a) ~ - cH 2 (t)(a-a min ) 2 . 



(72) 



Substituting this into Eq. (6), we can recast it as 

a + 3H(t)& + cH 2 (t)a = cH 2 (t) a min (t) , (73) 

where we see that there is now a source term on the right- 
hand-side (rhs). The solution to the above equation, 
therefore, is of the form 



0-(t) = (Thorn (t) + a pal (t) , 



(74) 



where the homogeneous solution a^ om {t) is given in 
Eq. (49) or Eq. (52) depending on the value of c. It 
is easy to see that the particular solution a pEiT (t) is of the 
form 



Cpar(t) = A a m i n (t) , 



(75) 



where, from Eq. (29), cr m i n (t) cx t "+ 2 and A is the con- 
stant 



A = 



l + 9 - 

c 



1 



n + 2 J Vn + 2 



1 



- 1 



(76) 



According to Eq. (71), the initial conditions, in this 
case, suggest that <7 cn d — c™^. Thus, a in Eq. (49) 
or Eq. (52) for the <7hom(t) is cr ~ (1 — A)a^^, because 
initially a(t ) = Atr™^. Therefore, for t » t cnd , the so- 
lution is 



a(t) ~ (1 - A)aT n — 



t 



4(^)( 1 -\^7^) 



+A a. 



end 



^end , 

t \ ~~ : 

^end , 



+ 



(77) 



where c = min{c, c x } and, in the oscillating case (i.e. 
when c > c x ), a(t) corresponds to the amplitude of the 
oscillations. 

It is interesting to note that the attractor constraint 
in Eq. (64) results in 



n <n r 



n+2 



< 



■5 (i+f) (i - V^W^j ■ (78) 



Therefore, we see that, if the attractor constraint is satis- 
fied, the particular solution decreases faster than the ho- 
mogeneous and, hence, it soon becomes negligible. If this 
is the case then there is no modification to the dynamics 
described in Sec. IV A other than the initial condition, 
which is the one given in Eq. (71). 

However, since the field is not in the higher order part 
of its potential, there is no compelling reason for the at- 
tractor constraint to be satisfied. Indeed, for sufficiently 
large n, the particular solution decreases less fast than 
the homogeneous. When this is so, the domination of 
the particular solution results in a(t) cx cr m i n . In other 
words, for sufficiently large n, a par takes over and keeps 
the field at a constant ratio A > 1 to er m ; n . 

The problem is that, in contrast to the homogeneous 
solution, the particular solution does not carry the mem- 
ory of the initial conditions and, therefore, does not pre- 
serve the perturbation spectrum of the field. This means 
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that the domination of the particular solution should be 
avoided, which imposes the constraint 

-l 



n + 2 < (n c + 2) 1- ,/l -mhJ 1 




(79) 



The above constraint is more relaxed than the attrac- 
tor constraint of Eq. (64) and becomes even more so the 
smaller c is. 

One final point to note before we conclude this section 
is that, from Eqs. (56) and (71), in the quasi-quadratic 
case with negative sign, after the end of inflation the 
curvature perturbation of the field is 
jj 

(a < C* < 1 . (80) 

Cmin on d 

We elaborate more on this point in Sec. VD3. 



V. CURVATON DOMINATION 

We now turn our attention to the evolution of the cur- 
vaton's density ratio over the background density Pa/p, 
which will determine when the curvaton will dominate 
the Universe, or, if it decays (or thermalizes) earlier, what 
is the value of r of Eq. (18). 

The history of the Universe may be quite complicated. 
There may be more than one periods of inflation, some 
of them lasting only a limited number of c-foldings (e.g. 
thermal inflation). There may also be more than one 
fields acting as curvatons, in the sense that they obtain 
a superhorizon spectrum of perturbations during infla- 
tion, which they manage, at some point, to impose onto 
the Universe. This is not improbable since the field con- 
tent of theories beyond the standard model is rich in 
potential curvatons. If this is the case, then the curva- 
ture perturbation of the Universe, that we observe today, 
corresponds e.g. to the curvaton which dominated the 
Universe last. (However, note that, if the curvaton that 
decays last only nearly dominates the Universe, there 
may be other significant contributions to the total cur- 
vature perturbation e.g. from other, earlier curvatons.) 
Being aware of these possible complications, we do not 
pretend to describe the evolution of the actual history 
of the Universe. Instead, we provide the tools to study 
this history. To demonstrate the use of those tools, we 
implement them in a minimalistic scenario, where there 
is one, long period of inflation and one curvaton field. 
In this minimalistic scenario, the inflationary period is 
followed by a matter era, dominated by the oscillating 
inflaton, and a radiation era, commencing after reheat- 
ing. We still ignore temperature corrections by assuming 
a negligible g postponing their study to Sec. VI. 

As mentioned in the previous section, the curvaton 
has a chance to dominate only if it decays (or thermal- 
izes) after the soft mass m dominates the effective mass 
m e fi = VV". Thus, we require 



where H m is the Hubble parameter when m ~ m c s, 
which is different in the quasi-quadratic case and the 
higher-order case. Below we treat both cases individ- 
ually again, as in the previous section. 



A. The quasi-quadratic case 

From Eqs. (50) and (53) and in view also of Eq. (14), 
we find 



P 



where 



K 




(82) 



(83) 



From the above and Eq. (48), it is evident that the expo- 
nent K depends on whether we are in the matter or the 
radiation era. Indeed, for the matter and radiation eras 
respectively, we have 



K 



MD 



and 



1 - y/l -min{l,16c/9} 
Krb = ^(l- Vl-min{l,4c}), 



which results in 



and 



Pa 




( #end / 


Pa 




P 


MD 




' P 


end 



Pa_ 




( — 


iCRD 

1 2z. 


p 


RD 


V r inf/ 


' p 



reh 



Note that, for c <C 1, we have 



K] 



MD 



- c and K RB 



c. 



(84) 
(85) 

(86) 
(87) 

(88) 



Thus, for c — > 0, the ratio p„l P remains constant as ex- 
pected according to the discussion after Eq. (54). 
Using Eqs. (10), (39) and (46), we find 



Pa 

P 



<~ c 
c„d V m P 



( Ccnd\ 2 > 1 / #cnd N 2 



mp 



and 



Pa 
P 



-Send / P 



(89) 



(90) 



end 



Hjn > r , 



(81) 



Now, in the quasi-quadratic case, it is easy to see that 
H m ~m/V^, (91) 
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Using this and Eqs. (86) and (87), we obtain 



P 



Cend 

~ C 

H~Hr- V m P 



x min < 1 



m/^/c 
H cn d 



r inf J 



• (92) 



When H{t) falls below H m , the quasi-quadratic term 
becomes dominated by the quadratic term correspond- 
ing to the soft mass m in Eq. (28). If c < 1 then there 
is a regime where something curious happens. Indeed, 
for m < H(t) < mj \fc, the scalar potential is dominated 
by the quadratic term, for which V" = m 2 <C H 2 , i.e. 
the field becomes overdamped and, consequently, freezes 
[70] . This does not happen in the case where the field was 
oscillating before H m was reached, because in this case 
c > Cx ~ 1. However, it happens when the field is follow- 
ing the scaling solution in Eq. (52), according to which, 
until H m was reached, the field was rolling steadily to- 
wards the origin. After H m , this roll is halted. The im- 
portant point to be stressed here is that, roughly when 
H m is reached, the field stops regardless of its value. That 
is the whole configuration of a freezes at the same mo- 
ment and, therefore, the perturbation spectrum is unaf- 
fected by the freezing. The reason is evident: the possible 
freezing has to do with the comparison of m and m e g(t) 
with each-other and with H(t). None of these quantities, 
in the quasi-quadratic case, is related to the actual value 
of g and, therefore, the freezing of the field configuration 
is a- independent. 

If the field configuration does freeze then, in the inter- 
val m < H < H m , the scaling of the curvaton's density 
ratio is only due to the decrease of p. Thus, from Eq. (13), 
in this interval, we have 



P 



(93) 



which is, for the first time, increasing. Is it possible for 
the curvaton to dominate (or nearly dominate) during 
this interval? Well, since p a j p oc p^ 1 oc H~ 2 , we find 



1 Pa_ 

C p 



(94) 



H~H„ 



Thus, after the whole interval of freezing has passed the 
amplification of the density ratio is only a factor of 1/c. 
Hence, at the end of the freezing interval, the curvaton's 
density ratio is 



Pa 

P 



)2 / / /—\Kmd { I r^RD-^MD 



where we have set (p a /p)m = {Pa/ 1 P)n~m- Notice now 
that, in the RHS of the above, all the factors are smaller 
than (or at most equal to) unity, while all the expo- 
nents are positive. This means that the density ratio 
of Eq. (95), even though somewhat increased compared 



to Eq. (92), is still expected to be much smaller than 
unity. It is hardly possible, therefore, that this ratio is 
able to satisfy the constraint of Eq. (25), let alone lead to 
domination. Therefore, curvaton domination will have to 
wait for the quadratic oscillations to begin. The initial 
density ratio for the curvaton at the onset of these final 
quadratic oscillations is given by Eq. (95). 

An important comment to be made here, regarding 
Eq. (95), is that the curvaton's density ratio (with re- 
spect to the background density at a given time) strongly 
depends on the value of c. Indeed, when c«l, Eq. (88) 
suggests that the two last factors on the RHS of Eq. (95) 
are negligible. On the other hand, if c ~ 1 then these 
factors are not negligible but, instead, they may strongly 
suppress the density ratio. Indeed, it is evident that the 
curvaton's density ratio, in this case, is suppressed at 
least by a factor of (m/i? en d), which can be very small. 
The reason is that, for c <C 1, the effective mass is small 
and the field is overdamped. On the other hand, if c 
is comparable to unity both the quasi-quadratic oscilla- 
tions and the scaling solution substantially diminish the 
density of the curvaton. In fact, it is possible that this 
suppression is so efficient that the final quadratic regime 
is unable to counteract it. Hence, we have shown that 
the quasi- quadratic evolution disfavours curvatons which 
admit substantial supergravity corrections to their effec- 
tive potential. This strengthens the case of the curvaton 
being a PNGB as discussed in Rcf. [41]. The above are 
more clearly demonstrated in the example analysed in 
Sec. VII. 



B. The higher— order case 

As discussed in Sec. IV B, in this case, the field un- 
freezes and begins oscillating when m e s ~ H. Suppose, 
at first, that, between the onset of the field's oscillation 
and the time when the soft mass m dominates m e ff, the 
dominant fluid component of the Universe remains the 
same so that w is unchanged. Then, using Eq. (66), it is 
easy to find that 



H„ 



( n+6 \( n c +2 n 
\n+2j\n c +e n 



±21 

+6) 



m . 



(96) 



where we used that m e s(H osc ) <~ H osc . Due to the at- 
tractor constraint in Eq. (64) we see that the exponent 
in the above is always positive. This means that 



H osc >m => H m < m . 



(97) 



The above shows that, the oscillation of the field will 
continue uninterrupted when the quadratic term takes 
over. This is because, for a pure quadratic potential, 
the field oscillations begin when H(t) ~ W" = m. Also, 
since H osc > H m (for the higher-order oscillations to take 
place), the condition H m < m is guaranteed if H osc < m. 
Therefore, in all cases H m < m and there is no interme- 
diate freezing of the field once it unfreezes after the end 
of inflation. 



13 



Now, employing Eqs. (63), (67) and (96), we obtain 



Pa 
P 

where we also used 



ra 



(l-™)(n + 6)-4 
re + 2 



^cnd 

ra P 
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(98) 



(99) 



which is evident in view of Eq. (10) and considering that 
Pa ~ ml s a 2 with m cS (H osc ) ~ H osc and p° sc = p c a nd . 

From the above, it is straightforward to find the cur- 
vaton's density ratio at the onset of the quadratic oscilla- 
tions in the case where this occurs before reheating (i.e. 
H m > Ti n f). Indeed, from Eq. (98), one gets 



Pa_ 

P 



MD 
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/ cr on d\ 
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corresponding to (c.f. Eq. (96)) 



H„ 



ra . 



(100) 



(101) 



Note that, for unfreezing before reheating, the attractor 
constraint demands n < 2. Similarly, we obtain the cur- 
vaton's density ratio at the onset of the quadratic oscilla- 
tions in the case when unfreezing occurs after reheating 
(i.e. H osc < r in f). This time, Eq. (98) gives 
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P 



H 



3(n + 2) 
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corresponding to (c.f. Eq. (96)) 
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ra 

Hqg ( 



3(n + 2) 



m . 



(102) 



(103) 



Note that, for unfreezing after reheating, the attractor 
constraint is relaxed to n < 6. 

What happens, however, if H m < Ti n f < H osc ? In this 
case, during the oscillations and before the quadratic 
term dominates, reheating is completed and the Universe 
switches from being matter dominated to being radiation 
dominated. We may obtain H m and (pa/p)m working in 
the same manner as above. 

Employing Eq. (66) both for the matter and radiation 
eras, we get 



ra 



H m 
rinf 



3/2 



H 



i + 2 
1 + 6 



H n 



(104) 



Solving for H m we obtain 

KtS) 

Hn 



Ho 



( ra \ 



1/3 



ra . 



(105) 



Since unfreezing occurs before reheating, the attractor 
constraint demands n < 2, which suggests that H m < m 
if H osc > Tti and ra < T ln {. Since H osc > H m , the con- 
dition H m < m is automatically satisfied if H osc < ra. 
Also, if ra > T; n f then again H m < ra is satisfied because 
H m < Tinf in the first place. Therefore H m < ra in all 
cases, which means that the oscillations continue unin- 
terrupted after the quadratic term dominates. 
Now, using Eq. (67), we get 
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(106) 
(107) 



Putting these together and also using Eqs. (99) and (105), 
we find 



Pa 
P 



ra \ 1 ^ rpT > 



Ho 



r \ 2/3 , N 2 

^ inf \ / ^end 



Ho 



rap 



. (108) 



Summing up, in the above we have obtained the 
value of H m [Eqs. (101), (103) and (105)] and of 
the ratio {pa/p) m [Eqs. (100), (102) and (108)] in 
the cases Tinf < H m < H OS c7 H m < _ff sc 5: Tinf and 
H m < Ti n f < H osc respectively. We have found that the 
attractor constraint ensures that H m < ra in all cases 
and, therefore, guarantees that the oscillations of the field 
continue uninterrupted after the quadratic term due to 
the soft mass m dominates the effective mass. 

Before moving on to the final quadratic oscilla- 
tions, we provide an estimate for -ff OS c- Indeed, since 



#osc ~ m 

gest that 



ff (a e „ d ) ~ <+7M 2 , Eqs. (37) and (39) sug- 



-ffosc > Vco-ffond • 



(109) 



C. The final quadratic oscillations 



When H(t) < H m , the scalar potential becomes dom- 
inated by the quadratic term due to the soft mass 



(110) 



According to Eq. (22), a scalar field coherently oscillat- 
ing in a quadratic potential corresponds to a collection of 
massive particles at rest (curvatons in our case) and be- 
haves like pressureless matter, i.e. p a oc a~ 3 . Therefore, 
from Eq. (13), one finds 



Pa 
P 



(111) 



The above shows that, in the quadratic regime the cur- 
vaton's density ratio is constant in the matter era and 
begins growing only after reheating. Thus, the curvaton 
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will have a chance to dominate (or nearly dominate) the 
Universe only if 



r < r inf . 

Using Eq. (Ill), one obtains 



Pa_ 

P 



~H 



1/2 



where 



Hfn = min{m, H m , T ini } 



(112) 



(113) 



(114) 



The above suggests that, if the curvaton decays (or ther- 
malizes) before domination, then, from Eq. (18), we have 



P 



(115) 



If, on the other hand, the curvaton does manage to dom- 
inate, then Eq. (113) gives 



Hdo 



(116) 



Hence, in the quasi-quadratic case, considering that 
Hm = min{m, Tinf} and in view of Eq. (95), the above 
gives 



#dom ~ min{m,r inf } 



. 4 / / i— \ 2K MD 
^end \ I m/ Vc 
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In the higher-order case there are three possibilities, as 
follows. The first possibility is that all the higher-order 
oscillations occur before reheating. Then, H m = r; n f and 
using Eqs. (100) and (116) we obtain 



Hdo 



2 / \ 4 

m \ I (Tend 



mp 



inf 



(118) 



The next possibility is that all the higher-order oscil- 
lations occur after reheating, in which case H m = H m . 
Then, substituting Eqs. (102) and (103) into Eq. (116), 
we find 



H 



dom 



H n 



°cnd 
77lp 



m . 



(119) 



Finally, the last possibility is that the higher-order os- 
cillations begin before reheating but they continue after- 
wards into the radiation era. Again we have H m = H m . 
Substituting Eqs. (108) and (105) into Eq. (116), we get 



Hdo 



H„ 



2 / \ 4 

fond 



mp 



rinf , 



(120) 



which is identical to the first possibility (c.f. Eq. (118)). 
Thus, in general, for the higher-order oscillations, we 
have found that 



Hdo 



H n 



( fond \ 4 
mi 



min{tf osc ,r inf }. (121) 



From the above, we see that, remarkably, in all cases 
Hdom is independent of n. 

What if the quadratic term was dominant right from 
the start? This is indeed possible if 



c < c < (--— ) « 1 . (122) 

\ -"end / 



Then, (p a /p) m = (p*/p) OS c, which is given by Eq. (99). 
Thus, in this case, we find 



-ffdc 



Cond 

mp 



min{m, r inf } . 



(123) 



Note that the above can be obtained by setting 
c (m/H cnd ) 2 in Eq. (117) and H osc -> m in Eq. (121). 



D. The value of q 

In this section, we take a closer look to the value of q, 
which relates the fractional perturbation of the curvaton 
at horizon crossing and at curvaton decay as defined in 
Eq. (4). We have already seen that, in all cases (c.f. 
Eqs. (56) and (69)), 



5a 
a 



5a 



end 



(124) 



Below, however, we will study this relation in more de- 
tail. For simplicity, we drop the bar from the oscillation 
amplitude of the field. 

During inflation and after the fast-roll regime, the cur- 
vaton is overdamped. As shown by Eq. (43), this means 
that the value of the curvaton perturbation is frozen. It is 
trivial to show, using Eq. (35), that a very similar relation 
is true for the value of cr* itself. This is especially true if 
the curvaton has reached the quantum regime. Thus, if 
c does not change significantly at the end of inflation, it 
is safe to consider that 



6a 
a 



5a 
a 



(125) 



end 



Also, since we have shown that, before decaying, the cur- 
vaton needs to be in the final quadratic oscillation regime, 
Eq. (21) suggests that 
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1. The quasi- quadratic case 

In the quasi-quadratic case, Eqs. (49) and (52) show 
clearly that 



a(t) =a x f{t/t ) , 



(127) 



where / is some function with /(l) = 1, which is indepen- 
dent of (To- It is clear from the above that, if we consider 
a perturbed value of the field, we will have 



6a 5a 

— (t) = = const. 

a cr 



(128) 



This is hardly surprising since, for potentials of the form 



V 



the equations of motion for the field and its 



perturbation (Eqs. (6) and (7) respectively) are iden- 
tical. Now it is crucial to note that the onset of the 
final quadratic oscillations occurs at the same time re- 
gardless of the value of the field. This is due to Eq. (91), 
which shows that H m is a-independent. Thus, choosing 
to = tend, Eq. (128) gives 



6 AH m ) = S - 
a a 



end 



Eq. (128) is obviously valid during the final quadratic os- 
cillations as well, carrying the perturbations intact until 
the time of curvaton decay, which again is cr-independent. 
Thus, 



5a_ 
a 



(130) 



dec 



As explained in Sec. VA, if c <C 1, it is possible that 
the field briefly freezes again before the onset of the final 
quadratic oscillations, but this is not expected to affect 
the amplitude of the perturbation spectrum. All in all, 
in view also of Eq. (125), we have found that, in the 
quasi-quadratic case 



which means that 



9=1 



2 So- 
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(131) 



(132) 



2. The higher-order case 

The case of higher-order term domination is a bit more 
complicated. Due to Eq. (65), one easily sees that, in 
this case too, Eqs. (127) and (128) are valid. However, 
because m c s ~ a n+2 /M n and in contrast to the quasi- 
quadratic case, the onset of both the higher-order and 
the final quadratic oscillations does depend on the value 
of a. Therefore, perturbing this value somewhat affects 
also the scaling of the field. Below we calculate what this 
means for q. 



Suppose we consider two causally disconnected regions 
of the Universe, characterized by two different values for 
the curvaton field: cr*- 1 ' and a^ 2 \ Suppose also that, 
in both these regions, the field lies in the higher-order 
regime. Then, in these two regions, the onset of the 
higher-order oscillations and the final quadratic oscilla- 
tions will occur at different times. Without loss of gen- 
erality, we assume that it is in region-(l) that the field 
begins its (Hubble damped) sinusoidal oscillations first. 
Using Eq. (65), we obtain the following for the ampli- 
tudes of the field at this moment: 



T (1) 
J m;l 



J2) 
J m;l 
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(133) 



Considering the fact that a oc H 2 /3(i+«>) an( j a i so 
#osc ~ ml s {a cnd ) oc cr™+, 2 , the above gives 



(129) 
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end 
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n±6\ 
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(134) 



Now, after the onset of the final quadratic oscillations 
for o-( x \ its density decreases as p a oc a 1 oc a~ 3 . On 
the other hand, still oscillates in the higher-order 
regime, until the final quadratic oscillations begin for 
too. Therefore, at this moment, we have 



°m;2 = O". 



T (2) 
7 m;2 



m;l 



7 (2) 
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(135) 



Let us assume now that, during the period of the 
higher-order oscillations, the equation of state of the Uni- 
verse remains unmodified. Then we can use Eq. (96) to 
relate H m with H osc . Using again that a oc H~ 2 / 3(1+w ~> 
and also that H 2 SC ~ ml s (a end ) oc tr™^ 2 , we find 
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m; 7 







(136) 



Now, after begins its final quadratic oscillations both 
fields oscillate sinusoidally (with Hubble damping) and, 
therefore, Eq. (128) is valid. Since the decay of the cur- 
vaton occurs independently of the value of <r, we see that 
(cr( 2 Ver (1) )dcc = (o-( 2 V<7 (1) )m;2- Using this and combin- 
ing Eqs. (134) and (136), we find that 
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dec 
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end 
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Therefore, if we take that a = a^ and a + 5a = a^ 
and also 8a <C a, we obtain 
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where the ellipsis denotes terms of higher power in 
(5a/a) cn d- Hence, in view also of Eq. (125), we have 
found that 



q= i(n + 6) 



1 



6 



n c + 2 



6 



n + 2 



(139) 



which, for matter and radiation eras respectively, be- 
comes 

2 — n 6 — n 

9md = — r~ and 3rd = -^T~ ■ ( 14 °) 



where c ~ 1 corresponds to the value after the end of in- 
nation. Thus, since a cnd (c* ) - a Q « <t££(c) < <C2(c), 
we see that the effective value of q may be substantially 
reduced if c changes significantly at the end of inflation in 
agreement with Eq. (80). Considering that cr ond (c*) w er» 
and also that 8a* ~ the above becomes 



<5a 



end 



CO 



From Eq. (44), wc have (5a /a)* 
in view of the above, gives 



(143) 

max{c*,c }, which, 



q c ff 



min < 1, 



Co 



(144) 



VI. THERMAL CORRECTIONS 



From Eq. (139), it is evident that, for n = n c , we 
have q = and all the perturbation spectrum is erased. 
From Eq. (137), it is clear that this is so to all orders in 
(8a/a) cn< \. Hence, we see that, the fact that the onset of 
the higher-order and the final quadratic oscillations is a- 
dependent suppresses strongly the perturbation spectrum 
if n = n c . This is why the border case n = n c for the 
attractor constraint in Eq. (64) is not acceptable, even 
though the field oscillates and does not converge towards 
the attractor solution. For n < n c , however, we see that 
q ~ 1 and so Eq. (124) is justified. 

In the case where the higher-order oscillations begin 
during matter domination but end after reheating, work- 
ing as above and using Eq. (105) instead of Eq. (96), one 
can show that 



9m D 



(n + 5)(2-n) 
3(n + 6) 



(141) 



Since, when unfreezing occurs during matter domination, 
only the case of the quartic term is allowed by the at- 
tractor constraint in Eq. (64), we see that, in this case, 
g=l/2 {q = 5/9} if the curvaton decays before {after} 
reheating. 



3. Effective q for significant change in c 

Before concluding this section, wc should point out the 
possibility that c does change significantly at the end of 
inflation. This will not have an effect if the field remains 
in the higher-order regime. However, if c ~ 1 after the 
end of inflation, this will probably relocate the field into 
the quasi-quadratic part of the potential. In this case, as 
discussed also in Sec. IV C, (5a/a) en< i may be substan- 
tially reduced. We will demonstrate this assuming that 
cr* has reached the quatum regime. 

In view of Eq. (71), one finds that Eq. (125) becomes 
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In this section, we introduce also the temperature cor- 
rection to the potential, as it appears in Eq. (28) . We will 
investigate the effect of this correction on the dynamics 
of the rolling and oscillating curvaton. The mass of the 
curvaton is taken now to include all the contributions 



V" 



m' 



±cH 2 (t)+m 2 B + g 2 T 2 



(145) 



where, in this section, with m e s we refer to the contribu- 
tion due to the higher-order terms: m 2 s <~ a n+2 /M 2 . 

During the matter era of the inflaton's oscillations, 
there exists a subdominant thermal bath due to the infla- 
ton's decay products. One would naively expect that any 
thermal bath during a matter era is to be rapidly diluted 
away since radiation density redshifts faster than matter, 
according to Eq. (13). However, due to the continuous 
decay of the inflaton field, the thermal bath receives an 
ever-growing contribution to its density, so that the con- 
tinuity equation, Eq. (12), does not apply and, therefore, 
neither does Eq. (13). The temperature of this thermal 
bath is [71] 



T < 



>r in f7?) 



1/4 



(146) 



According to the above, the radiation density during the 
matter era scales as p 7 oc T 4 oc H. Since this density is 
negligible compared to the one of the oscillating inflaton 
field, the loss of energy from the latter is also negligible so 
that, for the inflaton, the continuity equation does apply. 
Now, from Eq. (10), we have p oc H 2 , which means that 
there is a moment when the radiation density becomes 
comparable to the one of the inflaton. This is the moment 
of reheating and it occurs when -ff ron ~ r; n f. In view of 
Eq. (146), the reheat temperature is 
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\JmpT 



inf 
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After reheating, the Universe becomes radiation domi- 
nated with T 4 . In view of Eq. (10), therefore, we 
can write for the temperature of the Universe 
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Due to its interaction with the thermal bath, the scalar 
condensate is in danger of evaporating, that is of be- 
coming thermalized, if Tt > T a , i.e. if the decay rate 
of the field is smaller than the thermalization rate. If 
thermalization does occur, even if, strictly speaking, the 
curvatons have not decayed yet into other particles, they 
constitute a component of the thermal bath and their 
equation of state is that of radiation. Thus, after ther- 
malization the density ratio of the curvaton to the ra- 
diation background remains constant. As a result, r is 
decided not at the time of the curvaton decay but at the 
time of thermalization and this is why, so far, we have 
used r as given by Eq. (19) to determine the amplitude 
of the total curvature perturbation. 

Let us estimate Tt- The scattering cross-section for 
scalar particles is a ~ g 2 /B c 2 m , where E cm ~ \JTm a is 
the centre-of-mass energy. Thus, for the scattering rate 
r sc ~ ctT 3 , we find T sc ~ g 2 T 2 /m r7 . When the particles 
thcrmalize, m a ~ gT and so, for thermalization, we find 



gT. 



(149) 



Suppose, now, that the temperature corrections were in- 
deed dominating the effective mass of Eq. (145). Then 
V" ~ g 2 T 2 ~ T\. As usual, the field is able to oscillate 
only when H(t) drops enough so that V" ~ H 2 . How- 
ever, in this case, Eq. (149) suggests that, when the field 
is about to begin oscillating, H ~ Tt and, therefore, the 
field thermalizes and the condensate evaporates. Hence, 
we see that the field thermalizes before engaging into os- 
cillations, when the effective mass is dominated by the 
temperature corrections. 

As discussed in the previous section, the curvaton has 
a chance to dominate (or nearly dominate) the Universe 
only if it decays or thermalizes after the onset of the 
final quadratic oscillations. Thus, in order to protect the 
curvaton from thermalizing too early, we have to demand 
that 



gT m < m . 



(150) 



where T rn = T(H m ) is the temperature when VV" be- 
comes dominated by the soft mass to. Note that, for 
later times, gT < to and the thermalization rate is sup- 
pressed. 

The above condition, Eq. (150), suffices in order 
to avoid thermalization, only if, before H ~ H m , the 
temperature corrections have been always subdominant. 
This is ensured provided both cH 2 and m 2 s decrease 
with time faster than (gT) 2 . This is indeed the case as 
we show below. 

Consider first the quasi-quadratic case. Then, in view 
of Eq. (148), we find 
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which is indeed decreasing with the Universe expansion. 
Now, in the higher-order case, we have the following. 
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FIG. 2: Illustration (log-log plot) of the different scalings 
of the various contributions to V" assuming a given back- 
ground barotropic parameter w. The solid line corresponds 



to the higher-order case, where V" 



mitt 



a n+2 /M n . The 



dashed line corresponds to the quasi-quadratic case, where 
V" ~ cH 2 . The thin dot-dashed lines correspond to {gT) 2 
for different values of g growing as the thick dot-dashed arrow 
shows. It is depicted that the quasi-quadratic contribution to 
the effective mass decreases faster than the higher-order one 
and they both decrease faster than the temperature induced 
mass. From the figure, it is evident that, if (gT) 2 is subdom- 
inant when the soft mass-squared m 2 takes over, then this 
is enough to ensure that the temperature contribution to the 
effective mass is always subdominant, both before and after 
the moment when m 2 takes over. 



Eqs. (66) and (148), in view also of Eq. (63), give 
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The exponent in the RHS of the above is positive when 

n > —2 ( j+™ w J , which is always true for n > 0. Thus, 

the effective mass m e s also decreases faster than the gT 
and, hence, the condition in Eq. (150) ensures that the 
temperature corrections to the effective mass are always 
subdominant. This is illustrated in Fig. 2. 

The condition in Eq. (150) can be translated into an 
upper bound on the coupling g of the curvaton to the 
thermal bath. Indeed, from Eqs. (148) and (150), we 
find 



g 2 < 
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In the quasi-quadratic case, using Eq. (91), the above 
becomes 
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Now, using Eqs. (101), (103) and (105), in the higher- 
order case, we find 
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for r in f < H m < H osc , 

m / Hcy--\ 3(1+5) 
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(156) 



(157) 



for H m < H osc < Tinf . 

As discussed above, the constraint in Eq. (150) guar- 
antees that the temperature correction to the potential 
has negligible effect to the dynamics of the curvaton. In 
this case, T T becomes irrelevant and, therefore, T = T a . 
For strong coupling g, we expect T a ~ g 2 m. However, 
if g is very small then the curvaton may have to decay 
gravitationally, in which case its decay will be suppressed 
by mp. Thus, we can write 



r CT 



to max • 



m 
mp 



(158) 



VII. A CONCRETE EXAMPLE 

In this section, we will demonstrate our findings by em- 
ploying them in a specific model, under our minimalistic 
assumption of a single inflationary period and a single 
curvaton field. By assuming certain quantities, such as 
the inflationary energy scale, the decay rate of the in- 
flaton, the mass scale M and so on, we aim to obtain 
bounds on the curvaton model and its parameters, such 
as to or g. 

In particular, for our example, we choose the following 
values: 



M = m P , 



A„~l 



(159) 



and also 



V}{ 4 ~ 10 14 GeV , 
TO inf ~ 10 10 GeV , 



if* ~ 10 n GeV , 

r inf ~ lCT 6 GeV . (160) 



In the above choice, we took into account the COBE 
upper bound on the inflationary scale V^ 4 < 10 16 GeV 
necessary to render the inflaton's perturbations negli- 
gible [20]. Also, we considered that the mass of the 
inflaton m- m f should be small enough for slow roll, 
TO? nf <C H 2 , but not too small since, for most inflation 
models, r\ ~ (TOi n f/H*) 2 <~ 0.01. We chose mi n f to be 
of the so-called intermediate scale ~ ^/topto ow , where 
TO CW ~ 1 TeV is the electroweak scale. Finally, we as- 
sumed that the inflaton is a gauge singlet and, there- 



fore, decays only gravitationally so that T[ 



inf 



"inf 



/to 



The reheat temperature, for the above parameters, is 
T rch - 10 6 GcV (c.f. Eq. (147)), which is well above 
Tbbn but satisfies also comfortably the gravitino bound: 
T reh < 10 9 GeV [72]. 



We assume that i/end ~ if* and also that inflation 
lasts long enough for the curvaton to reach the quantum 
regime. Then Eq. (39) suggests 



(Tend ~ ctq ~ H*/ max{c„ c } . 



(161) 



This means that, in view of Eqs. (4), (18), (44), (56) and 
(69), one finds 



Co ~ rq max{c„ c } 



(162) 



where Co = 2 x 10 -5 . Using this, one can readily exclude 
the quartic case n — regardless of the choice of param- 
eters in Eq. (160). Indeed, reinstating the A„ in Eq. (40), 
one finds 



co 



;(Ao/3) 



1/3 



0.1 



(163) 



This means that, since c, « 1, we have c* < c . Thus, 
using Eq. (162) and employing also the WMAP constraint 
in Eq. (25) , it is easy to find that we require 



q < 10- 



(164) 



which is possible to realize only if q = q e s since, other- 
wise, q ~ 1. Thus, we need to have a change of sign for 
the quasi-quadratic term after the end of inflation. Now, 
using Eq. (144) in this case, we obtain 



q c ff 



(165) 



which can never be as small as 10~ 2 . Hence, the quartic 
case is excluded, unless A is extremely small or infla- 
tion ends well before the curvaton reaches the quantum 
regime. The remaining two possibilities are n = 2,4 (wc 
assume that symmetries in the potential allow only even 
powers). In the following, wc choose n = 4, because it 
offers the largest parameter space. Using this, Eqs. (40) 
and (160) give 



c ~ 10 



(166) 



We consider two cases. In the first case, we assume that 
the (approximate) symmetry which keeps c small during 
inflation persists after inflation ends, so that c = c, « 1. 
In the second case, we assume that c ~ 1 after the end 
of inflation, considering both the ±c cases. In all cases, 
we assume that c does not change significantly at the 
transition between the matter and radiation dominated 
epochs. 



A. Case c = c» ^ 1 

In this case, there is no change of sign for the quasi- 
quadratic term so we will drop the hat on a, since, when 
both the constraints of Eqs. (64) and (79) are satisfied, 
the curvaton's evolution does not depend on the sign in 
front of c. 
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Suppose, at first, that c < Co- Then the field lies into 
the higher-order part of the potential, where the attrac- 
tor constraint in Eq. (64) is applicable. For n > 0, we see 
that the curvaton must unfreeze after reheating, that is 



Tinf ^ Hose • 

Now, in view of Eq. (109), we find 

H osc ~V^ H *~ 10 9 GeV , 



(167) 



(168) 



where we used Eqs. (160) and (166). From Eqs. (160) 
and (168), it is clear that Eq. (167) is not satisfied and, 
therefore, the possibility of higher-order evolution (i.e. 
c < c ) is excluded. 

Let us assume now that cq < c <C 1, which results 
in quasi-quadratic evolution. Since q=l, in this case, 
Eq. (162) suggests that 



r ~ 10" 5 / c • 



(169) 



Using the WMAP constraint of Eq. (25), wc find 
that c < 10~ 3 . However, from Eq. (166), we have 
c > Co ~ 10 -4 , which allows only one possibility 



and, therefore, 



c~ 10" 



r ~ 10" 



(170) 



(171) 



which marginally satisfies Eq. (25). 

Now, since c<l, Eqs. (84) and (85) suggest that 
K M b ^ ^rd ^ (c.f. Eq. (88)). Then, using Eq. (95) in 
Eq. (115), we obtain 



mp 



rv 



1/2 



(172) 



Eq. (161) gives (T on d ~ 10 14 GeV. Substituting this into 
Eq. (172) and using also Eq. (160), we find the condition 



min{m, r inf } ~ 10 12 r, , 



(173) 



where we have also employed Eqs. (91) and (114). To 
study this, we need to use Eq. (158). Thus, wc consider 
the following two cases: 



• Small g: This case corresponds to 

,3 



9 < 



to 
mp 



m J 
mp 



(174) 



which means that the curvaton decays gravita- 
tionally. Using the above, the BBN constraint in 
Eq. (26) results in the bound 



> 10 4 GcV . 



(175) 



Hence, in view also of Eq. (160), we see that 
to S> Tinf . Consequently, Eq. (173) gives 



m ~ 10 6 GeV , 



(176) 



which satisfies the bound of Eq. (175). Inserting 
the above into Eq. (174), we obtain 



9 < 10- 12 , 



(177) 



which, in view of Eqs. (170) and (176), can be easily 
shown to satisfy the bound of Eq. (154). 



Large g: This case corresponds to 



> 



mp 



g m 



Using this, Eq. (173) is recast as 
minjl^j-iOV ■ 



(178) 



(179) 



Taking m < r in f , one finds g ~ 10 -6 , which can be 
easily shown, however, to violate the constraint of 
Eq. (154). Thus, we have 

m > r inf - 10~ 6 GeV . (180) 

In view of the above, Eq. (179) is now recast as 

g~ 10- 9 (m/GeV)- 1/2 . (181) 

Using this and Eq. (178), one can show that the 
BBN constraint in Eq. (26) is trivially satisfied. Fur- 
ther, Eq. (178) demands that 



m < 10 6 GeV 



(182) 



Finally, employing the constraint of Eq. (154), we 
find the bound 



to > 0.1 GeV . 



(183) 



Summing up, in this case, we have ended up with the 
following parameter space. For large g, we have found 



10 _1 GeV < to < 10 6 GeV, 
10~ 9 > g > 10~ 12 , 



(184) 



where the relation between to and g is given in Eq. (181). 
For small g, we have found 



to - 10 6 GeV, 
9 < 10 



-12 



(185) 



The above parameter space is shown in Fig. 3. Note that, 
for all the parameter space, the masslessness requirement 
to < H* - 10 n GeV is well satisfied. 



B. Case c ~ 1 3> c* 



In this case, Eqs. (162) and (166) yield 
rq < 0.1 . 



(186) 



The results depend on the sign of c. We treat the two 
cases separately. 
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FIG. 3: Illustration (log- log plot) of the parameter space 
(solid line) for a successful curvaton in the case when 
c = c* ~ 10 -3 , with model parameters given by Eq. (160). 
The upper bound g < 10~ 9 is determined by the requirement 
of Eq. (154), which ensures that (gT) 2 never dominates the 
effective mass. For g < 10 -12 , the curvaton decays predomi- 
nantly through gravitational couplings. 



1. The case of +c 



makes things much more difficult and, as a result, the 
parameter space is diminished. 

Now, employing Eqs. (160) and (190), one finds from 
Eq. (191) the condition 

m 2 - T CT x 7 2 10 38 GeV . (192) 

In view of the BBN constraint in Eq. (26), the above gives 
the bound 

m> 7 10 7 GeV. (193) 

Now suppose that g is small enough for Eq. (174) to 
be valid. Then, Eq. (192) gives m ~ 7~ 2 10~ 2 GeV, which 
strongly violates the bound in Eq. (193) for the whole 
range of 7. Thus, we have to consider the case when g is 
large and Eq. (178) is valid. Using this, Eq. (192) can be 
recast as 

g ~ 7 ~ 1 10- 19 (m/GeV) 1/2 . (194) 

In view of the above, the lower bound on g from Eq. (178) 
gives to < 7 _2 10 _2 GeV which again violates Eq. (193) 
for the whole range of 7. Hence, in this case there is no 
parameter space for a successful curvaton. 



In this case, q ~ 1 and, therefore, Eq. (186) suggests 
that 

10~ 2 < r < HT 1 , (187) 

where we also took into account the wmap constraint in 
Eq. (25). The above and Eq. (162) suggest that 

1 < 7 < 10 , (188) 

where we defined 

7 = max{l, c*/co} . (189) 

Using this and Eqs. (160) and (161), we obtain 

fTc„d-7 _1 10 15 GeV. (190) 

Now, when c ~ 1, it can be shown that, after infla- 
tion, the field finds itself in the quasi-quadratic regime. 
Eq. (187) shows that the curvaton has to decay before 
it dominates. For c > c x , Eqs. (84) and (85) show that 
Kmd = 1 and K^s — 1/2. Substituting this in Eq. (95), 
we find, from Eq. (115), that 




(191) 



By comparing the above with Eq. (172), we see that r 
is now smaller by a factor (m/i?*). This is because, in 
contrast to the case c -C 1 when the curvaton's density 
ratio Pal P remains almost constant after the end of infla- 
tion, in the case c ~ 1 this ratio is substantially reduced 
until the onset of the final quadratic oscillations. This 



2. The case of —c 

In this case, q = q e s as given in Eq. (144). Using this 
and Eq. (162), one, generically, obtains 

Co~rco(^p. (195) 

In view of Eq. (166) in our example, the above gives 

r ~ 1 (196) 

and, therefore, it is possible for the curvaton to dominate 
the Universe before its decay. Using Eqs. (144), (161) and 
(195), it is easy to show that, generically, 

CT min (c, H*) > (Tend <^ <7cff < 1 • (197) 

Thus, if <7 fF < 1, compared to the minimum after the 
end of inflation, the field appears to be at the origin. 
Therefore, the initial amplitude of the quasi-quadratic 
oscillations is 

(c,H,). (198) 

Using this in our example, Eqs. (29) and (160) suggest 

CTo„d ~ 10 15 GeV . (199) 

As before, for c > c x , Eqs. (84) and (85) show that 
Kmd — 1 and ifnrj = 1/2. Substituting this in Eq. (117), 
we obtain 

(200) 
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Inserting Eqs. (160) and (199) in the above, we find 

ffdom ~ 10- 40 (m/GcV) 2 GeV . (201) 

Since the curvaton decays after domination, we require 
that 

r CT < H dom . (202) 

Combining Eqs. (201) and (202) with the bbn constraint 
in Eq. (26), we find the bound 

to > 10 8 GeV . (203) 

Now suppose that g is small so that Eq. (174) is valid. 
Then, Eqs. (201) and (202) give to < 10~ 4 GeV, which 
strongly violates the bound in Eq. (203). Thus, we have 
to consider the case when g is large and Eq. (178) is valid 
instead. Using this, Eqs. (201) and (202) result in 

g < 10- 20 (TO/GeV) 1/2 . (204) 

It can be easily shown that, when the above is true the 
constraint of Eq. (154) is also satisfied. 

Now the lower bound on g coming from Eq. (178) is 

g > 10~ 18 (to/GcV) . (205) 

One can be easily check that this bound, in view of 
Eq. (178), suffices to satisfy also the BBN constraint in 
Eq. (26). 

The available parameter space is spanned by Eqs. (204) 
and (205) and it is a surface on the to — g plane, in con- 
trast to being a line, as shown in Fig. 3 for the previ- 
ous case. This is because, when the curvaton dominates, 
the necessary condition is the inequality in Eq. (202), 
whereas if it decays before domination the necessary con- 
dition is the equality in Eq. (115) (i.e. r has to have 
exactly the right value). 

Unfortunately, from Eqs. (204) and (205), we find that 
the parameter space exists only if to < 10~ 4 GeV, which 
strongly violates the bound in Eq. (203). Therefore, we 
have shown that there is no parameter space for a suc- 
cessful curvaton in the case when c ~ 1 regardless of the 
sign of c. This is due to the fact that the quasi-quadratic 
oscillations drastically reduce the energy density of the 
field. This reduction is impossible to be counteracted 
during the final quadratic oscillations. 

VIII. DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated and analysed the 
dynamics of the curvaton field. Apart from its soft mass 
to,, we have assumed that the field receives a contribution 
to its effective mass from supergravity corrections. This 
contribution is determined by the Hubble parameter and, 
hence, the potential acquires what we called the quasi- 
quadratic term. Another contribution to the effective 
mass may arise due to the presence of a quartic term or 



of non-renormalizable terms. Finally, we also considered 
possible thermal corrections to the potential, which arise 
due to the coupling of the curvaton to the thermal bath 
present after the end of inflation. 

Firstly, we discussed the behaviour of the field dur- 
ing inflation. We explained that the field is expected 
to engage, initially, in fast-roll, which soon, however, is 
terminated to be followed by slow-roll evolution. Even- 
tually, slow roll is also halted due to the action of field 
perturbations generated by quantum fluctuations. These 
perturbations effectively stabilize the coherent motion of 
the field for the remaining part of inflation. 

After inflation, we have studied the evolution of the 
curvaton's energy density in the cases when the effec- 
tive mass is dominated by the quasi-quadratic term or a 
higher-order one. In the case of the quasi-quadratic dom- 
ination, we have shown that there are two types of evolu- 
tion depending on how suppressed is the effective mass. 
For effective mass of order the Hubble parameter, we 
found that the field begins to oscillate, whereas, for sup- 
pressed mass, we have found a scaling solution, according 
to which the field rolls gently towards the minimum of its 
potential. Both solutions preserve the amplitude of the 
perturbation spectrum of the curvaton, obtained during 
inflation. 

In the case of the higher-order term domination, we 
have shown again that the field, even though frozen ini- 
tially, will, eventually, begin to oscillate (preserving also 
the amplitude of the perturbation spectrum) provided 
the order of the higher-order term is not too high. This 
is because, in the opposite case, after unfreezing, the field 
was found to follow an attractor solution, which loses all 
memory of initial conditions. Consequently, were the 
field to follow such an attractor, all the superhorizon 
spectrum of its perturbations would have been erased 
and the field could not act as a curvaton. 

A similar situation was also found in the quasi- 
quadratic case, when the effective mass was negative af- 
ter the end of inflation. In this case, if the order of the 
higher-order term (which decides the temporal position 
of the minimum of the potential) is too high, instead of 
rolling, the field was found to follow a particular solu- 
tion and remain at fixed, constant ratio with respect to 
the temporal minimum, losing thereby all memory of ini- 
tial conditions. The constraints obtained on the order of 
the principle higher-order term in both the above cases 
were comparable and their implementation was decisive 
for the evolution of the field. 

Generically, we have found the following behaviour for 
the curvaton after the end of inflation. Both in the quasi- 
quadratic case and the higher-order case, the energy den- 
sity of the field was found to be decreasing faster than the 
background density. Consequently, curvaton domination 
(or near domination) is possible only after the effective 
mass of the field becomes dominated by the soft mass. In 
all cases, when the attractor was avoided, we have found 
that the damping on the amplitude of the perturbation 
spectrum is negligible. Finally, we have shown that the 
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temperature corrections should not be allowed to domi- 
nate the effective mass of the field because if they do the 
field thermalizes before having a chance to dominate (or 
nearly dominate). This requirement introduces a strin- 
gent constraint on the coupling of the curvaton to the 
thermal bath. 

In all cases, we have calculated the density ratio r of 
the curvaton to the background density of the Universe 
at the time of the curvaton's decay. If the curvaton is to 
decay before dominating, this density ratio is necessary 
to determine the total curvature perturbation imposed by 
the curvaton field onto the Universe. We have also cal- 
culated when the curvaton dominates, if it does indeed 
so, in all cases considered. Furthermore, we calculated 
in detail the factor q relating the curvaton perturbation 
at horizon crossing with the one at its decay. Finally, 
we have employed our findings in a concrete model re- 
alization, for which we have shown that there is hardly 
any parameter space available for a successful curvaton 
if c ~ 1 after the end of inflation, regardless of the sign 
in front of c. This is due to the fact that, for large values 
of c, the density fraction of the curvaton is drastically 
reduced by either the quasi-quadratic oscillations or the 
scaling solution. As explained also in Sec. VA, we ex- 
pect this problem to be generic and, if so, this would 
favour models where c after inflation remains strongly 
suppressed, such as the PNGB curvaton models discussed 
in Ref. [41]. 

We should mention here that, due to the fact that our 
attempt to describing the curvaton evolution was aimed 
at the bulk of the parameter space, there are possibly 
special situations, not studied here, where the evolu- 
tion may be significantly modified. For example, if the 
quasi-quadratic term is dominant but negative after the 
end of inflation and the field is undergoing oscillations, 
there will be a moment, just before the soft mass taking 
over, when the quadratic and quasi-quadratic terms can- 
cel each-other. For a brief period, therefore, the field will 
be oscillating into a potential of higher order and, hence, 
will have an equation of state with a different (larger) ef- 
fective barotropic parameter as shown by Eq. (22). The 
outcome of all this will be a "glitch" to the scaling of the 
density ratio that may affect both the time of curvaton 
domination and/or the total curvature perturbation if the 
curvaton decays before domination. The importance of 
such an effect is possible to be studied only numerically. 

Another example, in the case of negative quasi- 
quadratic mass, is the effect of the local maximum at 
the origin, which may be felt if the initial condition of 
the field is tr en d ~ cr m i n , which is exactly what we expect 
in this case, if the sign of c changes at the end of infla- 
tion. Indeed, if the field engages into oscillations, there 



is a chance that, at some point, the oscillation amplitude 
marginally sends the field on top of the local maximum. 
The effect of such an event is the possible tachyonic am- 
plification of the perturbation spectrum. In an analogous 
situation, this tachyonic amplification (when the field, as 
it oscillates, is temporarily stabilized near a local max- 
imum of the potential) was thoroughly investigated by 
us in an earlier paper [6], where it was found that the 
amplification of the perturbations could be substantial. 
One could even imagine that such an amplification mech- 
anism may counteract the suppression of the amplitude 
of the perturbation spectrum due to attractor evolution. 

Furthermore, we can imagine a number of modifica- 
tions to our considerations arising due to a Universe his- 
tory more complicated than the one considered here. For 
example, one can insert a brief period of thermal inflation 
before or after curvaton domination. Another example 
would be to consider a period of 'kination', for which the 
background density has a stiff equation of state w ~ 1, 
which arises naturally in models of quintessential infla- 
tion [36, 37] or non-oscillatory inflationary models [34]. 
Indeed, the curvaton is a convenient mechanism to in- 
crease the effective reheating efficiency of such models 
[38] . We have not considered such possible complications 
in this paper mainly because there is an infinity of them. 
However, we tried to express our results in as much a 
model independent way as possible, so that they will be 
easy to implement in particular, more complicated sce- 
narios of the Universe history. 

All in all, we have investigated the evolution of the 
curvaton field during and after inflation. We have shown 
that a successful curvaton either oscillates or scales down 
to its minimum but always preserves the amplitude of the 
spectrum of its perturbations. We also showed that, in 
order to dominate (or nearly dominate), the curvaton has 
to decay after its effective mass becomes dominated by 
the soft mass term and the oscillation becomes (Hubble 
damped) sinusoidal. Furthermore, to avoid premature 
thermalization, the temperature corrections to the cur- 
vaton's potential should remain negligible throughout all 
its evolution. Finally, there is a part of the parameter 
space which results in destructive attractor evolution, 
which erases the perturbation spectrum and should be 
avoided. 
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